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ABSTRACT 

X-ray and optical observations of quadruply lensed quasars can provide a microarcsecond probe of the lensed 
quasar, corresponding to scale sizes of ~ 10 2 - 10 4 gravitational radii of the central black hole. This high 
angular resolution is achieved by taking advantage of microlensing by stars in the lensing galaxy. In this paper 
we utilize X-ray observations of ten lensed quasars recorded with the Chandra X-ray Observatory as well as 
corresponding optical data obtained with either the Hubble Space Telescope or ground-based optical telescopes. 
These are analyzed in a systematic and uniform way with emphasis on the flux-ratio anomalies that are found 
relative to the predictions of smooth lens models. A comparison of the flux ratio anomalies between the X-ray 
and optical bands allows us to conclude that the optical emission regions of the lensed quasars are typically 
larger than expected from basic thin disk models by factors of ^3-30. 

Subject headings: gravitational lensing — quasars: general 



1. INTRODUCTION 

In this paper we carry out a systematic, uniform study of ten 
quadruply gravitationally lensed quasars (hereafter "quads") 
for which one or more Chandra X-ray Observatory and op- 
tical images exist. We show how such observations can 
probe the lensed quasar on angular scales of a few micro- 
arcseconds. 

Simple models for the gravitational potentials of lensing 
galaxies are usually sufficient to reproduce the positions of 
4-image gravitationally lensed quasars. But these same mod- 
els - a monopole plus a quadrupole - fail to reproduce the 
optical fluxes of those images. Such "flux ratio anoma- 
lies" are thought to be the product of sma ll scale structure 
in the gravitational potentials of galaxies |Witt et aLlJ1995b 
| Mao & Schneiderll998tlC hiba 200 21 iMetcalf & Madaul200U 
balal & Kochanekl l2002t ISchechter & Wambsganssl 120021) . 
Of the two leading explanations, the more intriguing is that 
we are seeing the effects of dark matter condensations of sub- 
galactic mass. The more prosaic explanation (though exciting 
for very different reasons) is that the anomalies are largely the 
result of microlensing by stars in the intervening galaxy. 

A mass condensation can produce an anomaly only if the 
radius of its Einstein ring (a "circle of influence") is large 
compared to the emitting region of the lensed source. To set 
the scale for the discussion, th e Einstein radius o f a star in 
a typical lensing galaxy is ^3 (m/ Mq)(Gpc/ Dl) microarc- 
seconds, where Dl is the angular diameter distance of the lens 
and m is the mass of the star. If the optical continuum emis- 
s ion from the typical quasar originated from something like 
a Sha kura & Sunyaevl (1 19731) thin accretion disk, we estimate 
that the Einstein ring of a 0.7 M© star would be much big- 
ger than the optical emitting region. The stars in the lensing 
galaxy would then be expected to produce m icrolensing of the 
optical continuum (Mort onson et al.ll2005l) . As X-ray emis- 
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sion is generally thought to arise from an even smaller region 
than the optical continuum, the X-ray fluxes of quasar images 
ought likewise to exhibit such anomalies when observed with 
Chandra. 

This expectation for the X-rays is borne out by observa- 
tions. iMorgan et all d2001l) found that in the quadruple sys- 
tem RXJ09 11+0551 the flux ratio of the C image to the D 
image was 0.19 in the X-ray compared with 1 .28 in the op- 
tical. iBlackburne. Pooley. & R appaport d2006l) found that in 
the quadruple system RXJ1 131-1231 the A/B ratio was 0.18 
in the X-ray compared with 1 . 10 in the optical and model pre- 
dictions in the vicinity of 1 .70. More recent X-ray observa- 
tions show consid erable changes in the fl ux ratio s of the im- 
ages i n this lens (Koch anek et all f2006). And iPooley et all 
(2006) find that in the quadruple system PG 11 15+080 the 
AijA\ ratio is roughly 0.16 in the X-ray compared with 0.68 
in the optical and model predictions more nearly equal to 
unity. 

In each of the above cases, the sense of the anomaly is 
that the X-ray flux ratios are yet more anomalous (in the 
sense of disagreeing with the models) than the optical flux 
ratios. In previous work (cited above) we showed that this 
could happen only if the optical continuum emitting region 
were substantially larger than predicted for Shakura-Sunyaev 
disks. In a number of individual studies — th e three afore- 
menti oned cases, as well as MG J 04 14+ 05 34 (IChartas et al.1 
l2002l H 1413+ 117 (IChartas et alJl2004T) . and Q 2237+0305 
dDai et aD l2003) — microlensing has been invoked to explain 
some of the observed X-ray properties. In this paper we report 
the results of a systematic study of a larger sample of X-ray 
imaged quad lenses 5 . 

The above discussion illustrates how microlensing permits 
at least some resolution of a quasar on microarcsecond scales, 
two orders of magnitude better than VLBI. This corresponds 
to physical scales in the accretion disk of just a few thousand 
AU, or ^1000 gravitational radii for a~3x 10 8 M Q black hole 



5 T he systems are the quad ruple lenses HE 0230— 2130 IWisotzki et alj 
fT999l). MGJ0 414+0534 < Hew itt et all fl99l. RXJ0911+0551 
dBade et alj [1997I). SP SS J0924+0219 Hnada et all 12003d) PG 11 15+080 
j wevmann et alj 119801) . RXJ1131- 1231 ISluse et all 120031) . H 1413+117 
{MaeainetalJ 11983) . B 1422+231 JPatnaik et alj 119921). W FI J2033-4723 
Morgan et al . 2004), and Q 2237+0305 iHuchra e t all 19851) . 
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TABLE 1 
X-ray Fluxes and Flux Ratios 



Lensed Quasar 



Image Flux Ratios 



ObsID a 



Date 



Exp. (s) 



HS/HM HS/LM HM/LM LS/LM 



LM unabs. Fo.5-8keV c 
(10~ 14 erg cm" 2 s -1 ) 



HE 0230-2 130 



1642 2000 Oct 14.4 



14 764 



B/A 
0.44!! 



B/C 
0.701 



A/C 
1.6*5 



D/C 
4<j+ 008 



5.«ca 



MGJ04 14+0534 



417 
418 
421 

422 
1628 
3395 
3419 



2000 Jan 13.7 
2000 Apr 2.9 
2000 Aug 16.9 

2000 Nov 16.6 

2001 Feb 5.1 

2001 Nov 9.3 

2002 Jan 9.0 



6578 
7437 
7251 
7 504 
9020 
28413 
96 664 



A 2 / Ax 



0.50; 



,+0.1] 
0.10 



0.67*|* 
35+ 00 ' 
90 +0 09 

U -" W -0 OS 



A 2 /B A,/B C/B 



1.3 



+0.3 
-0.3 



1 G+0-4 
1.8-0.4 



9 £+0.4 



2 6 +05 



63 +0 11 

U.OJ_ 10 



65 +o n 



F B 

12% 

13!* 

15!* 

14!' 

16" 

n% 

14i 



RXJ091 1+0551 



419 
1629 



1999 Nov 2.7 

2000 Oct 29.8 



28 795 
9 826 



A/B A/D 



2.7!, 



4.6! 



B/D 



C/D 
0.35!°!? 



F D 
1.9 
1.8 



-0.6 
+0.5 
-0.4 



SDSSJ0924+0219 



5604 2005 Feb 24.0 17 944 

PG 11 15+080 

363 2000 Jun 2.8 26492 

1630 2000 Nov 3.3 9 826 



RXJ1131-1231 



4814 2004 Apr 12.2 10047 

H 1413+117 

930 2000 Apr 19.7 38 185 

5645 2005 Mar 30.1 88 863 



D/A 
1 4+ 007 



10+ 001 
u - lw -o.oi 



A/B 

1.8*5 



D/B 
0.45** 



62+ 013 

"• uz, -0.12 

1 2* 3 



A/B A/C 



22* 03 

-0.02 



A/C 



A/B 

•3 9+0.8 



A 2 /A x A 2 /C Ai/C 



B/C 
9 9+0.1 

B/C 

9+0.9 



C/B 

"•^-0.14 

B/C 
1 . l+oi 



D/C 
30 +0 03 

D/C 



Fb 

i O+0.8 



50!' 

T 7+0.7 



B 1422+231 



367 2000 Jun 1.6 28 429 
1631 2001 May 21.5 10651 
4939 2004 Dec 1.6 47 729 



B/A 



62 +0 ' 

u - u -0.08 



B/C A/C 
1.1*1 1.6*1 
1.4*: 



0.87!! 
0.95!! 



D/C 



u. i i_ 0] 

0.08*» 



1.7! 



37!? 
40+12 

33!' 



WFIJ2033-4723 



5603 2005 Mai' 10.1 



15 420 



A 2 /A, A 2 /B Ai/B C/B 



l.l! 



0.87! 



0.64! 



Fb 

3 5 +0 - 9 



Q 2237+0305 



431 
1632 



2000 Sep 6.7 

2001 Dec 8.8 



30 287 
9 538 



D/A 
17 +001 

ul '-0.01 



D/B 
n q 5+0.09 



A/B 
5.0!" 



C/B 
2 l* 2 



5 9 +1 3 
6.1!" 



il The observation identifier of the Chandra dataset. b HS = Highly magnified Saddle point; HM = Highly magnified 
Minimum; LS = Less magnified Saddle point; LM = Less magnified M inim um. See i|4] c The unabsorbed flux of the 
LM image is computed from the best fit power-law model described in j|2.1l 



at an angular diameter distance of 1 Gpc. 

In addition, as suggested above, lensed systems present 
unique opportunities to study not only the lensed object but 
also the lensing object. As we will discuss in a forthcoming 
paper, these same observations of X-ray flux ratio anomalies 
permit mea surements of the dark matter c ontent of the lensing 
galaxies (ISchechter & Wambsganssl[2004l) . 

In ^2] we discuss the analysis of the Chandra archival data 
for ten quads. In $3]we describe the properties of the selected 
group of lenses in the optical band. We also present a uni- 
form set of models for these lenses produced with the same 
software and for a common set of model assumptions. The 
flux ratio anomalies are compared between the X-ray and op- 
tical images in $4] In $5] we draw conclusions concerning 
the sizes of optically emitting regions in these ten sources. 



Finally, in $6] we summarize our findings. Throughout this 
paper we adopt a cosmology with £1a = 0.7, Q,m = 0.3, and 
// = 70kms- 1 Mpc" 1 . 

2. X-RAY OBSERVATIONS 

The data were downloaded from the Chandra archive, and 
reduction was performed using the CI AO 3.3 software pro- 
vided by the Chandra X-ray Center 6 . The data were repro- 
cessed using the CALDB 3.2.2 set of calibration files (gain 
maps, quantum efficiency, quantum efficiency uniformity, ef- 
fective area) including a new bad pixel list made with the 
acis_run_hotpix tool. The reprocessing was done with- 
out including the pixel randomization that is added during 
standard processing. This omission slightly improves the 

6 |http ://asc. harvard, edu] 
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FIG. 1 . — X-ray lightcurves of the high magnification pair of images in each quad. All observations in TablefT]were used to make these 0.5-8 keV lightcurves, 
with multiple observations of the same quad separated by hash marks. The time delay between the pair (from the models described in iJ5]and Table0 is given 
and shown as a thick horizontal bar in the cases where it will fit on the plots. 

point spread function. The data were filtered using the 
standard ASCA grades and excluding both bad pixels and 
software-flagged cosmic ray events. Intervals of strong back- 
ground flaring were searched for, and a few were found. In all 
cases, the flares were mild enough that removing the intervals 
would have decreased the signal-to-noise of the quasar im- 
ages since it would have removed substantially more source 
flux than background flux within the small extraction regions. 
Therefore, we did not remove any flaring intervals. The ob- 
servation IDs, dates of observation, and exposure times are 



given in Table Q] 

2. 1 . Determining X-ray Flux ratios 

For each observation, a sky image was produced in the 0.5- 
8 keV band with a sampling of 0."0246 per pixel. Because 
of the significant overlap of the lensed images (especially 
the close image pairs) in many cases, the intensities were 
determined by fitting to each sky image a two-dimensional 
model consisting of four Gaussian components plus a con- 
stant background. The background component was fixed to 
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TABLE 2 

Comparison of Gaussian Fitting to Aperture 
Extraction ofSDSS 1004+4112 



X-ray Image Flux Ratios 



Method 


A/B 


A/C 


B/C 


D/C 


Gaussian fit 
Extraction 


77+ 03 
77+ " 3 


96 +0 04 
n 03+004 

u -"-0.04 


1 "^-0.04 
i - z,u -0.05 


57 +0 03 
55^> m 



a value determined from a source-free region near the lens. 
The relative positions of the Gaussian components were fixed 
to the separations given in the CASTLES online database 7 , 
but the absolute position was allowed to vary. Each Gaussian 
was constrained to have the same full-width at half-maximum 
(FWHM), but th is value was allowed to float. The fits were 
performed using ICashl (11979|) statistics and the P owell mini- 
mization method in Sherpa (Freeman "etani200l . 

From the best fit 4-Gaussian model, the image flux ratios 
were calculated for the high magnification pair (saddle point 
and minimum; HS & HM, respectively) as well as for each 
image relative to the less magnified minimum (LM) image. 
The uncertainties on these ratios were determined with Sherpa 
via the projection command, which varies each ratio in 
turn along a grid of values while all other parameters are al- 
lowed to float to the new best-fit values. The results are given 
in Table Q] 

Because the Chandra point-spread function is only approxi- 
mately described by a Gaussian, we sought to test this method 
by utilizing a Chandra observation (ObsID 5794) of the large- 
separation quad SDSS 1004+41 12, for which all four images 
are well separated 8 . We extracted counts from the 90% encir- 
cled energy region of each ima ge, as determined by ACIS Ex- 
tract v3.94 (Br oos et"aT1 l2C)02). and formed a number of flux 
ratios. We also followed the above method of fitting Gaus- 
sians. The agreement in flux ratios is excellent (see Table|2]i. 

Finally, a spectrum of the LM image was extracted for each 
observation with ACIS Extract and fit in Sherpa via a sim- 
ple absorbed power l aw. The absorption consist ed of a fixed 
Galactic component (iDickey & Lockmanl ll990l) plus a vari- 
able component. This simple model provided an acceptable 
fit in all cases, and the additional absorption component was 
usually consistent with zero. The 0.5-8 keV flux of the unab- 
sorbed power law is given in Table Q] 

2.2. X-ray Variability 

As the numbers in Table[T]indicate, many of the flux ratios 
vary to some degree for the quads that have been observed 
multiple times. This may be due to varying degrees of mi- 
crolensing or to normal quasar variability combined with time 
delays among the images. In fact, variability plus a time delay 
could masquerade as a flux ratio anomaly. FigureQ]shows the 
X-ray lightcurves of the sum of the high magnification pair 
of images for each system, which are seen to be fairly con- 
stant in all systems; only small amplitude (factor of two or 
less) variability is observed, even in cases when the length of 
the observation exceeds the predicted time delay between the 
brightest images 9 . 

7 |http : / /www, cf a . harvard, edu/ castles /| 

8 SDSS 1004+4112 Hnada et al.1 12003H1 is lensed by a dark matter- 
dominated cluster of galaxies, and is unique enough that we did not include 
it in our sample of lenses. 

9 Predictions of time delays are subject to large uncertainties, especially 
for pairs of images with small separations. Indeed, the time delay for 



TABLE 3 
Optical Observations 



Lensed Quasar Optical magnitudes 



Obs. date 


HS 


HM 


LS 


LM 


Filter 


HE 0230-2130 


B 


A 


D 


C 




2002 Jul 29 


19.22 


19.02 


21.21 


19.59 


F814W a 


MGJ0414+0534 


A 2 


A] 


C 


B 




1994 Nov 08 


21.36 


20.43 


22.10 


21.24 


F814W a 


RXJ0911+0551 


A 


B 


C 


D 




2000 Mar 02 


18.38 


18.64 


19.36 


19.66 


F814W a 


SDSSJ0924+0219 


D 


A 


C 


B 




2003 Nov 19 


21.59 


18.69 


19.86 


19.52 


F814W b 


PG 1115+080 


M 


M 


B 


C 




2004 Feb 22 


15.86 


16.08 


17.68 


17.26 


i' c 


RXJ1131-1231 


A 


B 


D 


C 




2004 Apr 12 


17.43 


17.42 


19.72 


18.44 


I* 


H 1413+1 17 


A 


B 


D 


C 




1994 Dec 22 


Mil 


17.84 


18.15 


18.06 


F814W a 


B 1422+231 


B 


A 


D 


C 




1999 Feb 06 


15.85 


15.88 


19.68 


16.41 


F814W a 


WFIJ2033-4723 


A 2 


Ai 


C 


B 




2003 Aug 01 


19.14 


18.68 


19.41 


19.32 


i' e 


Q 2237+0305 


D 


A 


C 


B 




1999 Oct 20 


17.39 


15.92 


16.77 


17.21 


F814W a 



a See h ttp : / /www . cf a . harv ard. edu/castle s" IKeeton et al] 
|2006) C Relative magnitudes from Pooley et al. 12006); zeropoint 
from this work. d Sluse et U] J2006 ) e Morgan eTal] <2004l) 

For the rest of the analysis, we utilized the observation 
with the highest signal to noise ratio for the quads observed 
multiple times by Chandra. We chose not to average over 
multiple epochs in order to avoid averaging out variations 
due to changes in microlensing. We use ObsID 3419 for 
MGJ04 14+05 34, ObsID 419 for RX J091 1+0551, ObsID 363 
for PG 1 1 15+080, ObsID 5645 for H 1413+1 17, ObsID 4939 
for B 1422+231, and ObsID 431 for Q 2237+0305. 

3. OPTICAL IMAGES AND LENS MODELS 

We turned to the existing literature for optical data with 
which to compare our X-ray flux ratios. For each lens, we 

o 

used data near 8000A, either Sloan i , Cousins /, or HST 
F814W. An effort was made to choose the observations clos- 
est in time to the deepest Chandra observation. The dates of 
the observations, along with the optical bandpasses and the 
image magnitudes, may be found in Table [3] The images are 
arranged according to their magnifications and parity (see $4}. 
Under ideal circumstances, the X-ray and optical observations 
would have been made on the same day, in order to mini- 
mize systematic errors resulting from quasar variability and 
microlensing variability. But for most of these lenses, such 
contemporaneous observations have not been made. Three 
lenses have X-ray and optical observations separated by about 
6 to 10 years, three by 2 to 4 years, and four by 15 months or 

RXJ1 131— 1231 is n ow known to be greater than predicted by a factor of 
~ 12 iMorgan et al. 2006). This is mitigated by the fact that the optical and 
X-ray observations for this particular source were made on the same day. 
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FIG. 2. — Representation of the deviations from the models in X-rays and optical. Each of the three frames for a system is constructed by placing Gaussians 
at the relative positions taken from the CASTLES online database. The leftmost frame in each set has the intensity of each Gaussian set to unity. In the center 
frame, the intensities are set to the ratio of the optical flux (normalized by the optical rms as defined by eq.(l)) to the model flux (normalized by the model rms; 
see eq. (1)). The same is done for X-rays in the rightmost frame of each set. The same color scaling is applied to every frame. For aesthetic reasons, the FWHMs 
of the Gaussians are a constant fraction of the frame size; a 0"5 scale bar is shown at the bottom right of each "unity" frame, and this frame also gives the image 
names and image types (see {|4). 



less. One of these, RX J1 131-1231, w as observed in both 
bands on the same day (ISluse et al.ll2006l) . 

These delays between observations can add systematic un- 
certainty to the results. However, there are reasons to be- 
lieve that their effect is not a strong one. The general lack 
of strong quasar variability seen in X-rays (see 32.21 >. coupled 
with the limited success of campaigns to measure lens time 
delays (which rely on quasar variability), suggest that quasars 
do not often vary by the factors that would be required to ex- 
plain the flux ratio anomalies. The fact that RXJ1 131-1231 
has an extremely strong discrepancy between X-ray and opti- 
cal flux ratios despite simultaneous observations in both bands 
shows that time variability cannot fully explain the anomalous 
ratios. 

We used Keeton's (2001) Lensmodel software, vl.06, to 
model each of the ten lenses as a singular isothermal sphere 
(SIS) with an external shear. This model has seven free pa- 
rameters (lens strength, shear strength (7) and direction (0 7 ), 
and the positions of source and lens), making it overcon- 



strained by the ten input measurements (the positions of four 
images and the lensing galaxy). The position measurements 
were obtained from the online CASTLES database. The ob- 
served fluxes of the lens images were not used as constraints. 

The models fit the image positions fairly well in all cases 
except that of HE 0230-2 130, where the position of the D 
image is significantly altered by a second galaxy. Since this 
lens has an obvious strong perturbation from a companion 
lens galaxy, we added a second mass component to the model. 
Allowing its position and strength to vary, and using its mea- 
sured position as a constraint, gave us eleven free parameters 
and twelve constraints. We found that a steeper projected pro- 
file than isothermal was required for this second mass compo- 
nent, so we modeled it as a circular power-law profile with an 
index of -1 .3. This model allows a much better fit to the data, 
and predicts an Einstein ring radius of 0."79 for the main lens- 
ing galaxy and 0"42 for the perturber, and an external shear 
of 0.10 in a direction 60.1° west of north. 

Parameters for the remaining lenses may be found in Ta- 



6 



Pooley et al. 



TABLE 4 
Lens Models 



Magnification b 

E 7 </> 7 11 HS HM LS LM zi z s 



HE 0230-2 130 c 










-11.80 


11.50 


-2.32 


6.22 


0.52 


2.162 


MGJ0414+0534 


1" 


20 


0.13 


77.1° 


-20.72 


19.07 


-1.68 


5.36 


0.96 


2.64 


RXJ091 1+0551 


1" 


11 


0.32 


1.1° 


-4.41 


8.10 


-3.23 


1.77 


0.77 


2.80 


SDSSJ0924+0219 


0" 


88 


0.04 


84.6° 


-23.19 


26.78 


-12.57 


10.98 


0.39 


1.524 


PG 11 15+080 


1" 


15 


0.12 


65.0° 


-13.37 


14.54 


-3.02 


3.88 


0.31 


1.72 


RXJ1131-1231 


1" 


78 


0.12 


-73.3° 


-23.72 


13.93 


-1.58 


13.40 


0.295 


0.658 


H 1413+117 


0" 


61 


0.11 


21.8° 


-5.17 


5.46 


-3.32 


5.05 




2.55 


B 1422+231 


0" 


78 


0.27 


-54.6° 


-12.04 


8.86 


-0.35 


5.69 


0.34 


3.62 


WFIJ2033-4723 


1" 


12 


0.15 


26.3° 


-6.61 


7.69 


-2.20 


3.17 


0.66 


1.66 


Q 2237+0305 


0" 


88 


0.07 


67.1° 


-9.81 


9.21 


-5.32 


7.93 


0.04 


1.69 



a Measured in degrees East of North. b Negative magnifications signify saddle point 
images. c HE0230-2130 has a unique mass model with an extra companion galaxy. See text 
for lens parameters. 



ble |4] The predicted magnifications may be expected to 
vary with different ch oices of lens models at the 10% level 
dMetcalf & Zhaoll2002h . 

Our model, in which the quadrupole term of the gravi- 
tational potential arises from an external tide, gives larger 
magnifications (and therefore smaller bolometric luminosi- 
ties) than would a model in which the quadrupole is due to the 
flattening of the lens galaxy. Holder & Schechterl (120031) have 
argued that the high ratio of quadruply lensed quasars to dou- 
bly lensed quasars can be explained if most of the quadrupole 
is tidal in origin. 

4. COMPARISON OF ANOMALOUS FLUX RATIOS: 
X-RAY VS. OPTICAL 

Figure|2]provides a visual guide to the optical-to-model and 
X-ray-to-model flux ratios of each quad. It shows represen- 
tations of each system using two-dimensional Gaussians, the 
positions of which come from the CASTLES database. As 
a point of reference, the leftmost frame for each quad shows 
Gaussians of unit amplitude. The center frame represents the 
optical-to-model ratio of the images, normalized by each rms 
(described below). The amplitude A,- of image i is given by 

A i = ^f^ (1) 

I Mi I A ms | n | 

where i = 1,2,3,4, F opt)I is the (linear) optical flux of image 
i, and //, is the image magnification from Table [4] The right 
frame gives a similar representation for the X-rays. The rms 
of the optical (and X-ray) observations is first computed as 




However, because the rms can be dominated by one highly 
anomalous image, we remove the largest deviator and then re- 
compute the rms. The largest deviator is defined as the image 
i with the maximum value of | log 10 (A,'. opt ) + log 10 (A,, x - ray )|. 
This new rms is then used in eq. (1) to compute the ampli- 
tudes, the values of which are given in Table [5] 

In every case save one, the most anomalous image was the 
highly magnified saddle point image. This is not surprising, 
since ISchechter & Wambsganssl d2002l) have shown that mi- 
crolensing is likely to affect high-magnification saddle points 



most strongly. In order to give the lenses a uniform treatment, 
we have classified the four images in each lens according to 
their magnifications and the local morphology of the travel- 
time surface. Henceforth in this paper, "HS" will designate 
the highly magnified saddle-point and "HM" the highly mag- 
nified minimum. Likewise, "LS" will designate the less mag- 
nified saddle-point and "LM" the less magnified minimum. 

In this work, we are most interested in the optical-to-model 
and X-ray-to-model ratios of the HS/HM flux ratio. The com- 
parison between optical and X-ray ratios is shown for each 
quad in Figure [3] The first panel shows the observed HS/HM 
ratio relative to the model HS/HM ratio, and the second and 
third panels show how each of HS and HM compare to the 
less magnified minimum image (LM). In almost all cases, the 
HS/HM ratio is more extreme in X-rays than in the optical; 
when the observed ratio is greater then the model ratio, the 
X-ray ratio is greater than the optical, and, when the observed 
ratio is less than the model ratio, the X-ray ratio is less than 
the optical. The second and third panels show whether the 
discrepancy with the model comes from the HS or the HM 
image (or a combination of the two). In general, the LM im- 
age is much less susceptible to microlensi ng than either the 
HS or HM image dKochanek & Dalalll2004b . 

The group statistics for the flux ratio anomalies presented 
in Figure [3] and Table [5] are summarized in Figure |U The 
error bars represent the ±rms spread in the logarithm of 
the flux ratios (normalized by the smooth model values) be- 
tween various image pairs for our quasar sample. The black 
outer bars result from including all 10 quasars; the heavy 
blue bars result when we exclude the systems Q 2237+0305 
and SDSSJ0924+0219. Q 2237+0305 is excluded because 
the uniquely small redshift of its lensing galaxy causes the 
projected microlens Einstein radius to be bigger than any 
region of the source, while SDSS J0924+0219 might also 
be excluded because the source size is thought to be so 
small that eve n its broad line region is partially microlensed 
dKeeton et al.ll2006l) . 

It may be seen from the blue bars in this figure that the 
ratios of the HS to HM images deviate more (from their ex- 
pected values) in the X-ray band than in the optical band by a 
factor of —2.4. The discrepancy is somewhat smaller for the 
HS/LM ratios at a factor of -1.7. The HM/LM and LS/LM 
ratios are not as anomalous in either band, but the X-ray ra- 
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FIG . 3 . — Comparison of X-ray (blue X ) and optical (red O ) ratios to lens model ratios for select image pairs for each lensed quasar. The leftmost frame shows 
the ratio of the highly magnified saddle point (HS) to the highly magnified minimum (HM), while the center and rightmost frame show the ratio of each of these, 
respectively, to the less magnified minimum (LM). The ratios for the X-ray are based on the observation with the highest signal to noise, and those for the optical 
are based on the observation closest in time to the chosen X-ray data. The light blue X's show the variation in the X-ray ratios for quads observed multiple times 
by Chandra. 



tios still have a wider range than do the optical ratios. It is 
on the larger anomalies in the X-ray band for the HS/HM and 
HS/LM ratios, as compared to those for the optical band, that 
we base our quantitative analysis of the size of the optically 
emitting regions of the accretion disks in the next section. 

5. SIZES OF QUASAR EMISSION REGIONS 

For the purpose of interpreting our results, we adopt the 
working hypothesis that the anomalous flux ratios presented 
in this paper are the result of microlensing. Microlensing by 
stars in the lensing galaxy can account for the observed flux 
ratio anomalies, but only if the source is small compared to 
the Einstein radii of the microlensing stars. Figure [3] shows 
dramatic evidence for microlensing in the X-ray band for at 
least 7 of the 10 lensing systems in our study. In general, 
the optical emission of these same systems, while still being 
microlensed, has less extreme flux ratio anomalies than in the 
X-ray band by a factor of ^2 (see Figure|4]and the discussion 
above). Since the X-rays are expected to be emitted very near 
to the black hole, the condition for microlensing is easy to 
meet - the source should indeed be quite small compared to 
the Einstein radius of the microlensing stars. By contrast, the 
markedly lower degree of microlensing in the optical band 
implies that the size of the optical emission region in many of 
these sources is roughly comparable to the size of the stellar 
microlens Einstein radius. 

Many authors have studied the effect of source size on the 
microlensing of quasars by intervening galaxies. Typically 
the results are presented as plots of microlensing light curves 
(e.g., Wambsganss & Paczynski 1991) rather than rms fluc- 
tuations in the logarithm of the flux. There are no analytic 
techniques for estimating rms fluctuations, so one must simu- 
late the microlensing process. 

Ideally we would run point source simulations for each of 
the 40 images in our sample, taking into account the theoret- 
ical magnification (which in turn depends upon two indepen- 
dent parameters, a convergence and a shear) and the fraction 
of baryonic matter. Each simulation would produce a magni- 
fication map, which might then be convolved with sources of 
different sizes, producing magnification histograms. 

Such an effort lies beyond the scope of the present paper, 

10 While the rms microlensing fluctuations in flux are formally divergent, 
rms fluctuations in the logarithm of the flux are not (e.g., Witt et al. 1995). 



but we can draw upon such simulations that have been carried 
out. In particular we use the work of Mortonson et al. (2005) 
who studied in detail the effect of source size on minima and 
saddlepoints with magnifications of +6 and -6 respectively, 
assuming that the convergence (a dimensionless surface den- 
sity) is due entirely to equal mass stars and taking the con- 
vergence to be equal to the shear, as would be the case for an 
unperturbed isothermal sphere. The magnifications observed 
for our highly magnified minima and saddlepoints are larger 
than this, typically by a factor of two, while our less magni- 
fied images are typically fainter than this by a factor of two. 
Moreover there is reason to think that the stellar component 
comprises only a fraction - somewhere between 1/10 and 1/2 
- of the mass surface density. In the absence of a complete set 
of simulations we take those of Mortonson et al. as represen- 
tative. 

They find that, independent of the detailed radial profile 
of the source, the rms logarithmic fluctuations depend only 
upon the ratio of the half-light radius of the source to the Ein- 
stein radius, rin/r-Em- The rms logarithmic fluctuations de- 
crease from their maximum value at ri/i/r-Em = to one half 
that value at rx^l^m ~ 1/3. Since our optical fluctuations 
are roughly one half the amplitude of the X-ray fluctuations 
(which we take to arise from a region of negligible extent) we 
infer that the line-of-sight projected size of the optical region 
is roughly 1/3 the Einstein radius of the stars. 

To estimate a rough size for the expected region of the opti- 
cal emission from quasar accretion disks, we adop t a generic 
thin disk model (see, e.g., Sha kura & Sunvaevlll973l) . In such 
a model the gravitational energy release is redistributed via 
internal viscous stresses in such a way that, independent of 
the detailed nature of the origin of the viscosity, the rate of 
energy release per unit area of the disk at radius, r, is: 



(3) 



where M, M, and rg are the black-hole mass, accretion rate, 
and the inner radius of the accretion disk, respectively. Note 
that in this formulation neither special nor general relativistic 
effects are included, except implicitly via the location of ro. 
In our context, such relativistic effects are unimportant in the 
case of a Schwarzschild black hole. Relativistic corrections, 
including those for accretion disks around Kerr black holes 
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TABLE 5 

Flux-to-model Ratios Normalized by RMS 



HS/HM 



Optical 



X-ray 



n- 





Band 


HS 


HM 


LS 


LM 


HE 0230-2130 


Optical 


0.76 


0.94 


0.62 


1.02 




X-ray 


0.46 


1.06 


1.50 


1.24 


MGJ0414+0534 


Optical 


0.43 


1.10 


2.68 


1.86 




X-ray 


0.58 


1.03 


2.33 


1.72 


RXJ0911+0551 


Optical 


1.49 


0.64 


0.82 


1.14 




X-ray 


1.76 


0.36 


0.24 


1.29 


SDSSJ0924+0219 


Optical 


0.09 


1.10 


0.80 


1.25 




X-ray 


0.20 


1.20 


0.34 


0.92 


PG 1115+080 


Optical 


1.02 


1.15 


1.04 


1.19 




X-ray 


0.20 


1.14 


1.48 


1.10 


RXJ1131-1231 


Optical 


1.03 


1.76 


1.87 


0.72 




X-ray 


0.10 


1.72 


2.08 


0.82 


H1413+117 


Optical 


1.12 


0.99 


1.23 


0.88 




X-ray 


1.23 


0.69 


0.96 


0.87 


B 1422+231 


Optical 


0.76 


1.00 


0.77 


0.96 




X-ray 


0.49 


1.22 


1.78 


1.10 


WFIJ2033-4723 


Optical 


0.71 


0.93 


1.65 


1.25 




X-ray 


0.84 


0.65 


1.66 


1.81 


Q 2237+0305 


Optical 


0.33 


1.36 


1.08 


0.48 




X-ray 


0.22 


1.41 


1.00 


0.33 



NOTE. — The rms values were computed from the three least 
anomalous images in each quad. See 0for details. 



(e.g., iPage & Thornd[l974l) are only likely to exacerbate the 
difficulties with understanding the size of the optical emission 
regions discussed below. 

In the context of the thin-disk model around a 
Schwarzschild black hole, the fractional luminosity that 
emerges within a radial distance r is 



f L «r) = 3r / 1 



r~ l dr = 1 - 



.3/2 



'fa (r/r ) 

(4) 

The complement of this quantity, {I- fi(< r)}, i.e., the frac- 
tion of the luminosity released at radii > r, is plotted in Fig- 
ure [5] Here we have labeled the axis in physical units start- 
ing at ro = 6GM/c 2 = 2.5 x 10 14 cm, i.e., the last stable or- 
bit about a Schwarzschild black hole of 3 x 10 8 M Q , an illus- 
trative quasar mass. We also show curves for other possible 
black hole masses. For black holes with appreciable angular 
momentum, the value of ro moves progressively inward, and 
radii at which equal fractions of the luminosity are emitted do 
likewise. 

Also overplotted on Figure [5] are nine arrows, one for each 
of our sources with known redshifts, marking the physical size 
of the Einstein radius of a 0.7M Q star in the lensing galaxy as 
projected back onto the lensed quasar. What we see is that 
the arrows are virtually all located at radii where only a tiny 
fraction of the quasar luminosity can emerge from the disk - 
at least for our fiducial black-hole mass of 3 x 1O 8 M . These 
fractional luminosity values are typically < 2% for sizes com- 



HS/LM H 1 — 

HM/LM H- 
LS/LM H- 



-0.5 0.5 -0.5 0.5 

RMS Deviations in log(Observed Ratio/Model Ratio) 

FIG. 4. — The rms of the flux ratio anomalies in the optical vs. X-ray 
(see Fig.f3]and Table 5 for the flux ratios of individual sources). The black 
outer bars result from including all 10 quasars in our sample; the heavy blue 
bars result from excluding Q 2237+0305 and SDSS J0924+0219; the red bars 
result when we exclude only Q 2237+0305. 

parable to the backprojected stellar Einstein radii. Only for 
black-hole masses > 3 x 10 9 M Q does a significant fraction 
of the luminosity (i.e., ^20%) originate from radial distances 
comparable in si ze to the Einstein ra dius. However, even then, 
as we showed in lPoolev~et al. (2006) and demonstrate below, 
much of this radiation should be emitted at wavelengths well 
beyond the optical or near IR. Given that the optical radiation 
(e.g., 0.4-1.5 pm) typically compri ses a substantial f raction 
of quasar luminosities, e.g., ^15% (lElvis et al.lll994l) . it ap- 
pears difficult for the optical emission to be released from a 
thin disk at radii that are sufficiently large to allow for the par- 
tial suppression of microlensing — as observed. We further 
quantify this conclusion below. 

Figure [5] and eq. (4) imply effective upper limits to the size 
of thin accretion disks in the optical by evaluating the bolo- 
metric luminosity emitted within a radial distance r of the cen- 
tral black hole. We now proceed to compute more quantita- 
tively how large the accretion disk is expected to appear for a 
fixed waveband, e.g., V, R, I. Based on the relativistic invari- 
ant 7„/V 3 ,we find the following expression for the half-light 
radius, r^n, of a thin accretion disk in a waveband centered at 
v (in the Earth's frame): 



r i/2 [ hu(l+z)/kT(r)_ 



rn/2 

Jin 



1 rdr 



J™ [ e Mi+z)/«>)_i] l rdr 



(5) 



where ro is the location of the inner edge of the accretion disk, 
and T(r) is the local temperature of the accretion disk, which 
in the Shakura & Sunyaev (1973) model, is 



T(r) = 



3GM m M 



87rerr 3 



1/4 



(6) 



In this simple picture, calculation of the half-light radius re- 
quires knowledge of three parameters: M bh, M, and rp. We 
use primarily the optical-based method of Kaspi et al.l d2000l) 
(discussed below) to estimate the bolometric luminosity of 
each of the 10 sources in our sample. We also utilize the X-ray 
luminosity, coupled with a bolometric correction factor (also 
discussed below) to provide a sanity check on the Kaspi et 
al. approach. We further assume that all of the quasars are op- 
erating at the sa me fraction, Ie — 1 /4 , of their respective Ed- 
dington limits (Koll meier et al.ll2006l) . We show below from 
a simple scaling argument, that our final results for Tin are 
relatively insensitive to this choice. Finally, we assume that 
the radiation efficiency (rest mass to radiant energy conver- 
sion efficie ncy, 77) of all the quasar s in our sample is r\ = 0.15 
(see, e.g., lYu & Tremaindl2002l) . For this choice of effi- 
ciency, the dimensionless black-hole spin parameter would 
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Disk radius r (cm) 

FIG. 5. — Fraction of luminosity (fij emitted beyond radius r in a geo- 
metrically thin accretion disk for a variety of black hole masses. The arrows 
indicate the physical sizes of the Einstein radii of 0.7 Mq stars in each of the 
nine lensing galaxies of known redshift projected back onto the the lensed 
quasar. 

be a = 0.88 and the innermost stab le orbit would be located 
at r ~ 2.5R g = 2.5GM Bli /c 2 (e.g., lBardeenlll970t) . However, 
in our simple non-relativistic disk model, we can only fix ro, 
and accept whatever the non-relativistic energy release is. For 
ro = 2.5 R g this turns out to yield an equivalent r\ = 0.2, which 
is sufficiently close to the Kerr value to provide the desired 
accuracy in computing r\ n. 

We summarize the computed and inferred properities of 
our quasar sample in Table [6] The second c olumn gives 
the bo lometric luminosity as calculated from the iKaspi et all 
(2000) prescription. In this approach, Lboi is taken to be 
9[AFx] 51 oo471y/£. To estimate the 5100A flux in the rest 
frame of the quasar, we used the flux measured in the closest 
available broadband filter, usually the HST NICMOS F160W 
band, and e xtrapolated using an as sumed power-law spectrum 
F\ - A" 1 7 (iKollmeier et all 120061) . The third column in Ta- 
ble [6] gives an independent estimate of the bolometric lumi- 
nosity for each quasar based on the measured X-ray lumi- 
nosity and a bolometric correction f actor of 20, as in ferred 
from the composite AGN spectrum of lElvis et al.l(ll994l) . The 
fourth column provides the black-hole mass inferred from the 
bolometric luminosity (in column 2) divided by the Edding- 
ton fraction, /g = 1/4, which then yields Lndd, and thence 
Mbh- It should be noted that since Kol lmeier et al . (2005 ) 
derive f E = \/A using the prescription of IKaspi et alj (2000), 
and since we follow suit, the masses we derive are indepen- 
dent of the dimensionless factor (of 9) in Kaspi's prescrip- 
tion. The error bars on the mass represent the uncertainties 
inferred from the ± rms (logarithmic) spread in the bolomet- 
ric luminos i ties o btained via three different estimates: (i) the 
Kaspi et al. (2000) method, (ii) the X-ray luminosity, and (iii) 
(in 7 of the 9 cases) the mass estimate directly provided by 
Peng et al. (1999; based on a "virial" method involving broad- 
line widths and sizes of broad-line regions). The values of the 
half-light radius, r\n, computed with eqs. (5) and (6) for the 
I band are given in the fifth column. In the sixth column are 
Einstein ring radii of typical O.7M microlenses, projected 
onto the plane of the source. Finally, in the last column we 
give the logarithm of the ratio of the half-light radius to the 
microlens Einstein radius. 

The results of our thin-disk estimates for the ratio r\jijr-E m 
(last column of Table |6]l are plotted in Figure|6] The central 
heavy point within each error bar is based on the black-hole 
mass given in column 4 of Table|6] The error bars on the ratio 
r i/i/ r Em are propagated from the black hole mass uncertain- 
ties given in Table [6] The vertical line at r^/tin = 1/3 is 
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FIG. 6. — Ratio of thin-disk half-light radii to typical (0.7 Mq) microlens 
Einstein radii for the ten sample lenses. Estimates of the minimum value 
required to significantly attenuate microlensing variability are 0.1 (red) and 
0.33 (blue) - see text for an explanation. 

our estimate of the ratio required to suppress microlensing in 
the optical band by the (logarithmic) factor of ^2 discussed in 
§|4] The vertical line at r\ n/rsm =1/10 represents a more con- 
servative lower limit on rxn/rEm that might plausibly still be 
consistent with the suppressed microlensing in the optical. An 
inspection of Figure|6]shows that 7 of the 9 systems (for which 
r \l2/ r Vm could be calculated) lie below the limit of 1/10, and 
therefore the disk size in the optical that we calculate appears 
to be too small to explain the reduced microlensing. Only 
one of the systems, B 1422+231, has a ratio of rin/r^m that 
slightly exceeds the 1/3 value which we think is reasonable 
to account for the reduced microlensing in the optical for this 
particular source. 

Arguably the largest contribution to the uncertainty in the 
calculation of r^n arises from the errors in estimating the 
bolometric luminosities of the quasars in our sample. We be- 
lieve we can make a fairly robust estimate of the uncertainty 
in r\n — due to the uncertainties in Lboi — by inspection of 
Figure[6] As discussed above, the plotted error bars are de- 
rived from the logarithmic rms scatter among the three (two) 
different and independent methods we have employed to infer 
Lboi (see above discussion) for 7 (2) of the sources. As Fig- 
ure|6]indicates, the uncertainties in r\ n range between factors 
of 1.2 and 3.2, with an average value of a factor of 1.7. We 
take this to be a fairly reliable estimate of the uncertainty in 
our values for r\n due to errors in estimating L DO i. 

In our calculations leading to the set of values for ri/2/ r Em 
we assumed values for two key parameters of the quasars: (i) 
the radiative efficiency 77, and (ii) the fraction /e = Lboi/Lndd- 
Based on a simple scaling argument, we can show how our 
results for r\n depend on 77 and fy. Equation (5) provides the 
exact definition of r x > 2 that we use. However, if we use the ex- 
pression for T(r) in eq. (6) to find the radius where T /(l+z) 
equals hv/k, where v is the center of the observation band, 
this is to a good approximation proportional to r x n. If we fur- 
ther neglect the factor (1 - y/ro/r) in eq. (6), we find a handy 
scaling relation for rm: 



r 1/2 oc (M BH M 



1/3 



(7) 



If we consider the bolometric luminosity to be a measured 
quantity for each system, then Mbh °c A>oi//b and M oc 
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TABLE 6 
Quasar Properties 



Quasar i*oi,opt a £boi.x b log Mbh c ri/j d r 1/2 d stellar r Ein c log ryz/rmn 

(lO^ergs -1 ) (10 45 ergs-') (M ) (10 15 cm) (R g ) (1() 15 cm) 



HE 0230-2130 


2.9 


6.3 


7.95 ±0.24 


0.93 


70 


43 


-1.66±0.16 


MGJ04 14+0534 


36 


28 


9.04 ±0.17 


3.8 


23 


31 


-0.91 ±0.11 


RXJ0911+0551 


13 


13 


8.60 ±0.1 8 


1.9 


32 


35 


-1.26±0.12 


SDSSJ0924+0219 


0.6 


0.3 


7.27 ±0.56 


0.42 


152 


48 


-2.06 ±0.37 


PG 1 1 15+080 


11 


6.6 


8.53 ±0.37 


2.5 


50 


55 


-1.35 ±0.25 


RXJ1131-1231 


0.80 


1.3 


7.39±0.19 


0.84 


230 


38 


-1.65 ±0.13 


H 1413+117 


56 


6.5 


9.24 ±0.51 


5.4 








B 1422+231 


250 


135 


9.89±0.18 


13 


11 


47 


-0.55 ±0.12 


WFIJ2033-4723 


5.7 


3.8 


8.24±0.12 


1.6 


62 


36 


-1.35 ±0.08 


Q 2237+0305 


32 


2.7 


8.99 ±0.76 


5.5 


38 


150 


-1.43 ±0.51 



a Bolometric luminosities computed using L^„[ = 9[\F\]5iQQ4ird^. Computed from HM, LS, and LM images, corrected 
for magnification (Kaspi et al. 2000). b Approximate bolometric luminosities derived from the X-ray (0.5-8 keV) lumi- 
nosities (computed from LM image) with a bolometric correction factor of 20 (see ^5). c Calculated from the bolometric 
luminosities in column 2. See d rj n is computed according to eq. (5) for the I band. c Einstein radius of a O.7M0 star, 
projected back to the lensed quasar, in units of 10 15 cm. 



Lboi/v- Combining these, we can see how r\n depends on 
the assumed parameters and r\: 

rin^ifEVT 1 " , (8) 

which is a fairly weak dependence, and not likely to lead to 
uncertainties in r^n of more than an additional factor of ~2. 

6. SUMMARY AND CONCLUSIONS 

We have presented a study of ten quadruply gravitation- 
ally lensed quasars for which high spatial resolution X-ray 
and optical data are available, paying particular attention to 
the differences between the observed flux ratios of the high 
magnification pairs of images (i.e., HS/HM) and the pre- 
dicted flux ratios from smooth lensing models. The Chandra 
data were analyzed in a uniform and systematic manner, and 
the X-ray flux ratios were determined via two-dimensional 
Gaussian fits. The optical fluxes and image positions were 
found in the existing literature, with the bulk coming from the 
CASTLES project. We also modeled each lensing system as 
a singular isothermal sphere with external shear (except for 
HE 0230-2 130, where a second mass component was neces- 
sary), and these simple models fit the image positions quite 
well. 

As illustrated in Figures |2]-|4] almost all systems show ev- 
idence for an anomaly in the ratio of high-magnification sad- 
dle point and minimum images (HS/HM) as compared to the 
smooth model prediction. In the systems which show a pro- 
nounced anomaly, the X-rays are generally seen to be more 
anomalous than the optical. 

For a number of reasons, we believe that the anomalous 
flux ratios, and the differences between these ratios in the X- 
ray and optical ban ds, are best explained by microlensing . 
In previous work dBlackburne. Poolev. & R appaport l2006t 
iPoolev et al.l l2006) we have shown that extinction in the vis- 
ible band and absorption of soft X-rays cannot provide the 
explanation. Second, we show in this study (as well as pre- 
vious work) that temporal variability intrinsic to the source, 
in conjunction with lens time delays, also cannot, in most 
cases, explain the observed anomalies. Third, since images 
in both the X-ray and optical bands exhibit these flux ratio 
anomalies, but to differing degrees, no smooth lens model 



can reproduce these anomalies. Finally, we find that in the 
preponderance of systems, it is the highly magnified sad- 
dle point image (HS) whose flux is anomalous. This is 
in agreement with microle nsing magnification distributions 
(ISchechter & Wambsganssll2002l) . Since there is no reason 
for the HS location to systematically produce larger optical 
extinctions or X-ray absorptions, this is another argument 
against differential extinction/absorption being the cause of 
the flux ratio anomalies. 

Under the hypothesis that the anomalies are produced via 
microlensing by stars (of typical mass 0.7 M Q ) in the lens- 
ing galaxy, the implication is that the optical emitting region, 
which suffers rms (logarithmic) microlensing variations only 
half as big as those of the X-ray region, must have a typical 
size ~l/3 of the Einstein radius of the microlensing stars (see 
discussion in Likewise, the X-ray emitting region, being 
more severely microlensed, must be substantially smaller than 
this. 

In the context of a thin accretion disk around a black hole, 
the X-ray requirement is easily satisfied, as this emission 
likely arises from the inner parts of the disk. However, the 
optical emission poses something of a problem. It is generally 
thought to arise from a region not much larger than the X-ray 
region, but this is in conflict with the observed microlensing 
results which require larger optical emitting regions by factors 
of ^ 3 - 30 (see Figure|6]l than are commonly accepted. 

Therefore, we are left with a conundrum. Either there is 
a mechanism to transport the optical radiation to larger radii 
(and which does not affect the X-rays), or there is a miss- 
ing piece of the puzzle. Regardless, we have demonstrated 
how the X-ray and optical observations can provide a micro- 
arcsecond probe of the lensed quasars, and thereby yield po- 
tentially important results. 

From the work in this paper and the above discussion we 
draw three summary conclusions: 

• microlensing is the primary cause of the flux ratio anoma- 
lies. 

• the optical emitting regions in the quasars involved in this 
study have sizes of ~l/3 of a stellar Einstein radius, i.e., ~ a 
microacrsecond, corresponding to ^1000 AU. 

• millilensing (e.g., by dark matter haloes) is ruled out as an 
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explanation of the flux ratio anomalies by virtue of the above 
conclusion since this implies that both the X-ray and opti- 
cal emission regions are small compared to the milliarcsecond 
scale, and should therefore be lensed by the same amount. 

Finally, as mentioned earlier in the paper, these same flux 
ratio anomalies can be used to provide valuable information 
on the ratio of stellar matter to dark matter in the lens galaxy in 
the vic inity of its Einstein radius. ISchechter & Wambsganssl 
(2004) used optical flux ratio anomalies for a sample of eleven 
quads to derive a projected stellar/dark mass ratio at the typ- 
ical impact parameter of a quasar image. They first assumed 
that the optical emission region was small compared to the 
Einstein rings of the stars in the lensing galaxy. The re- 
sult was very heavily influenced by the inclusion (or exclu- 
sion) of the system with the most extreme flux ratio anomaly, 
SDSS J0924+0219. They found less discordant results if they 
instead assumed that half the optical light came from a point- 
like source and half came from a more extended source. In 
this case the results indicated a stellar to dark matter ratio of 
~ 1/4 (at the typical impact parameter of a quasar image), with 
an uncertainty of a factor of about 2. But allowing for a frac- 
tion of the light to come from a more extended region adds a 
second parameter to the problem, making determination of the 



stellar/dark matter ratio more uncertain. If, as we have argued 
here, the X-ray emission comes from a region substantially 
smaller than the optical emission region, the use of X-ray flux 
ratio anomalies in the analysis of Schechter & Wambsganss 
would eliminate that second parameter and more uniquely de- 
termine the ratio of stellar matter to dark matter. The same 
study would give a better idea of the emission region size re- 
quired to attenuate microlensing variations. Such an analysis 
is the subject of a forthcoming paper. 
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